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difficult  to  determine,  117  records  from  single-charge  events  appeared  (with 
simple  corrections)  to  give  credible  peak  velocities;  368  gave  times-of- 
arrival.  Inconsistencies  among  the  records  reduced  these  numbers  to  73  and 
29 2,  respectively. 

The  arrival  times  give  a  wavespeed  of  4.479±0.021  km/s,  which  differs  insig¬ 
nificantly  from  the  Cowboy  speed.  In  addition,  with  the  yield,  W,  in  kilotons 
peak  velocities  decay  with  r  according  to  the  expression 

0.24(r/Wa)  •  ^m/s,  over  the  scaled  range-interval  0.38<  r/W3  <11.3 

km/kt^.  Thus,  peak  velocity  is  nearly  twice  that  observed  in  Cowboy  where 
the  two  range-intervals  overlap,  but  it  decays  at  a  slightly  slower  rate. 

Most  important:  The^Cowboy  Trails  data  show  clearly  that  peak  velocity  decays 
inelastically  (pCr  )  even  at  7  times  the  largest  previous  range  of  observed 
radial  motion  from  tamped  explosions  in  dome-salt. 

Further  study  of  the  velocity  gauge  to  determine  the  cause  of  baseline 
problems  and  correct  for  them,  plus  the  application  of  digital  processing 
procedures,  should  allow  recovery  of  useful  data  from  the  bulk  of  the  flawed 
records.  Independent  measurement  of  displacement  (attempted  in  one  of  the 
events  considered  here)  would  provide  a  direct  test  of  correction  accuracy. 
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SECTION  1 


INTRODUCTION 


Careful  analysis  of  ground-motion  data  from  the  5.3-kt 
Salmon  nuclear  event1  and  the  series  of  tamped  chemical  explosive 
events  conducted  during  the  Cowboy  program  has  raised  basic 
questions  about  the  seismic -wave  sources  produced  by  explosions 
in  salt,  and  about  the  properties  of  the  explosive  -  pelletized 
TNT  ("Pelletol")  -  used  in  the  Cowboy  experiments3.  Specifically 

1.  Decay  of  peak  velocity  and  peak  displacement  with  slant  range 
show  that  salt  behaved  in  a  markedly  inelastic  way  out  to  the 
largest  scaled  ranges  instrumented  in  Salmon  (0.487  km/kt^)  and 
Cowboy  (5.744  km/kt3).  At  what  range  and  amplitude  does  defor¬ 
mation  become  elastic  (i.e.,  far-field  peaks  ocR-1)? 

2.  Wave  propagation  can  be  linear  without  being  elastic.  The 
Salmon  and  Cowboy  measurements  don’t  permit  linear  inelastic  be¬ 
havior  to  be  identified,  but  seismic  theory  assures  us  that  at 
"low  enough"  amplitudes,  geo-materials  deform  in  linear  fashion. 
At  what  slant  range  does  deformation  become  linear  in  dome-salt, 
and  what  does  the  explosively  generated  linear  source  look  like? 

3.  If  Pelletol  releases  the  energy  usually  attributed  to  TNT, 
then  the  curves  of  peak  velocity  and  displacement  vs .  range  for 
Cowboy  lie  below  those  for  Salmon.  What  is  the  correct  energy 
release  for  Pelletol  charges,  tamped  in  dome  salt? 

4.  With  what  accuracy  does  simple  scaling  apply  to  motion  from 
tamped  explosions  in  salt? 

The  four-Phase  Cowboy  Trails  program  was  undertaken  under 
Defense  Advanced  Research  Projects  Agency  (DARPA)  sponsorship  to 


answer  these  questions.  Phases  I  and  II  were  conducted  to  deter¬ 
mine  the  energy  release  of  tamped  chemical  explosive  charges  in 

4  5 

dome-salt,  by. measuring  the  resulting  cavity  volumes  ’  -  and  also 

to  evaluate  simple  scaling  rules  in  terms  of  those  volume^.  Phases 
III  and  IV  of  the  Cowboy  Trails  program  were  designed  i)  to  find 
where  the  response  of  dome-salt  to  a  burst  in  it  becomes  sensibly 
linear,  and  what  the  motion  is  like  there,  ii)  to  locate  any  trans¬ 
ition  from  inelastic  to  elastic  propagation  of  outgoing  signals, 
and  iii)  to  determine  the  accuracy  of  simple  scaling  rules.  Appen¬ 
dix  A  contains  a  detailed  account  of  the  program  plan. 

Peak  material  velocity  in  fields  driven  by  single  spherical 
charges  has  been  taken  so  far  as  the  criterion  for  elastic  propa¬ 
gation:  At  far  ranges  (i.e.,  where  distance  to  shot -point  is  » 
wavetrain- length  from  first  arrival  to  peak  velocity)  peak  velocity 
would  decay  inversely  with  slant  range  in  a  homogeneous,  isotropic 
elastic  medium.  As  for  evaluating  the  rules  of  simple  scaling, 
it  suffices  to  compare  simply-scaled  properties  of  fields  produced 
by  bursts  of  as  widely  different  yield  as  possible.  In  this  pro¬ 
gram,  measurements  were  made  of  the  sizes  and  shapes  of  cavities 
left  by  chemical  charges  weighing  as  little  as  20-30  gm;  by  con¬ 
trast,  cavity  measurements  are  available  from  bursts  as  large  in 
yield  as  Salmon.  As  regards  scaling  of  motion  and  its  Fourier  com¬ 
ponents,  extensive  ground-motion  data  were  acquired  in  this  program 
for  charges  weighing  ~890  N;  again,  comparable  data  extend  to  yieldc 
as  high  as  Salmon's,  so  that  pertinent  yields  span  almost  a  factor 
of  forty  in  the  scale  of  time  and  distance.  In  the  matter  of  source 
linearity,  no  comparisons  with  earlier  events  are  needed  (and  there 
are  no  relevant  in-situ  measurements  from  earlier  events  in  dome 
salt).  Rather,  the  measure  of  linearity  is  the  accuracy  of  super¬ 
position.  Thus,  the  program  included  events  in  which  two  spherical 
charges  at  various  separation-distances  were  simultaneously  detona¬ 
ted,  so  that  the  resulting  superposed  motions  (which  were  measured) 
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could  be  compared  with  single-charge  sums. 

In  view  of  Cowboy  data,  the  effects  of  principal  interest 
in  this  program  were  sought  at  greater  scaled  ranges  than  those  of 
prior  free-field  measurement  in  salt.  Hence,  comparisons  with 
Salmon,  especially,  could  be  made  only  at  the  smaller  of  the  ranges 
observed  in  these  shots.  Other  complications  stem  from  a)  the  use 
of  a  dome  not  previously  employed  for  ground-motion  experiments, 
b)  use  of  spherical  charges  rather  than  Cowboy-like-cylinders, 
and  c)  the  exclusion  of  Pelletol  from  the  smallest  bursts  of  the 
program  (due  to  its  large  grains).  The  logical  scheme  for  evalu¬ 
ating  these  non-ideal  factors  one  by  one  is  rather  intricate; 
Appendix  A  goes  into  its  details. 

SECTION  2 

THE  EXPERIMENTS 

The  planned  program  of  ground  motion  experiments^  at  the 
Morton  Salt  Company's  Kleer  Mine  near  Grand  Saline,  Texas,  inclu¬ 
ded  i)  an  event  with  a  single  spherical  charge  at  the  center  of 
the  instrumentation  array  (Event  III,  comprising  Phase  III  of  the 
program) ,  and  ii)  a  series  of  three  charge-pair  events  (Events 
IV-1,  IV-2,  and  IV-3,  comprising  Phase  IV)  with  charges  located 
symmetrically  about  the  main  line  of  ground-motion  instruments. 

Instrumentation  for  the  Phase  III  and  IV  events  consisted 
of  a)  68  velocity  gauges  and  three  accelerometers  aligned  with  the 
r,  8,  and  ^-directions  of  a  spherical  coordinate  system  with 
origin  at  the  Event  III  charge-center,  and  b)  17  accelerometers 
and  7  two-axis  displacement  gauges  aligned  with  the  R,  Z,  cp-direc- 
tions  of  a  cylindrical  coordinate  system  with  symmetry  axis 
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passing  vertically  through  the  Event  III  charge  center.  The 
spherically  aligned  gauges  were  deployed  at  eleven  slant  ranges 
along  four  lines  (denoted  A,  B,  C  and  D)  corresponding  to  rays 
from  the  center  of  a  regular  tetrahedron  (the  Event  III  charge) 
through  its  vertices,  as  shown  in  Figure  1.  The  velocity  gauges, 
which  formed  almost  all  of  the  primary  gauge-array,  were  electron¬ 
ically  integrating,  Be 11 -and -Howe 11  piezoelectric  accelerometers 
of  three  types:  Model  4-155-0122  for  the  closest-in-stations , 

Model  4-155-0111  at  one  intermediate  range,  and  Model  4-155-0129 
for  the  larger  slant  ranges.  Sunstrand  Model  305B  servo  acceler¬ 
ometers  replaced  the  tangential  (ep)  velocity  gauges  at  stations 
11,  12,  and  13  on  the  A-line  (hence  denoted  All,  etc.).  Two  types 
of  accelerometer  were  used  in  the  back-up  ground-motion  array 
along  line  A:  Columbia  Research  Laboratories  Model  904-PC  at  ' 
stations  3,  4  and  5,  and  Columbia  Research  Laboratories  Model  3029 
adjacent  to  stations  5,  6,  7,  9,  11,  12  and  13.  The  two-axis 
displacement  gauges  (experimental  freo-fall  gauges  being  developed 
by  PAI^)  were  installed  at  stations  A3,  A4  and  A5 ,  B3  and  B5 ,  C3, 

and  D5  for  use  only  in  Event  III.  Under  separate  contract,  S-Cubed 

8 

fielded  surface  seismometers  for  some  of  the  events  and  Stanford 

9 

Research  Institute  fielded  Flatpak  gauges  to  measure  stresses  in 
Event  III  near  the  charge. 

Delays  in  reaming  spherical  cavities  for  charge  placement, 
problems  in  demonstrating  the  safety  of  drilling  deeper  than 
177  m  below  the  working-level  floor,  and  misfires  in  the  charge- 
pair  events  led  to  the  measurement  of  ground  motion  for  the 
ten  events  (in  firing  sequence)  in  Table.  1. 

Charge  locations  for  the  ground  motion  events  are  shown  in 
Figure  2.  Nominal  charge  weight  was  890N  (200  lb)  but  actual 
charge  weights  varied  from  672N  to  850N  (151  to  191  lbs).^  Four 
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II 


a)  ELEVATION:  A-D  PLANE 


Figure  1.  Ground  Motion  Ins trumentation  Stations  for  the  Cowboy  Trails  Events. 
Except  as  noted  canisters  contained  three  orthogonally  oriented  (r,0,!p)  velocity 
gauges.  Accelerometers  were  substituted  for  cp-velocity  gauges  at  Stations  All, 
A12  and  A13.  On  lines  B,  C  and  D  (which  contain  only  odd  numbered  stations) 
canisters  at  Stations  5,  9  and  13  contained  only  r  and  g  velocity  gauges.  Radial 
(R)  accelerometers  were  installed  along  line  a  at  Stations  3,  4,  5,  6  and  7,  an 
R,Z-pair  at  Station  9  and  orthogonal  triads  (R,Z,cp)  at  Stations  11,  12  and  13. 

No  gauges  were  installed  at  B13,  D3,  or  D13. 
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Figure  2.  Plan  View  of  the  Ground-Zero  Region  Showing  Charge  Locations  (★)  for 
the  Cowboy  Trails  Ground  Motion  Events.  Open  circles  indicate  a  number  of 
instrument  stations  to  facilitate  comparison  with  the  instrumentation  layout  of 
Figure  1.  Close-in  gauge  locations  (A3,  AA,  A5,  B3,  C3,  and  C5)  have  been 
omitted  for  clarity 


of  the  single-charge  events  (IV-1.0,  IV-1.1,  IV-2.0  and  IV-2.1) 
were  the  result  of  misfires  during  attempts  to  simultaneously 
detonate  symmetric  charge  pairs.  One  charge,  IV -lb  failed  to 
detonate  in  events  IV-1.0  and  IV-1.1,  the  first  two  attempts  to 
fire  the  dual  charge-pair  ever.t  IV-1.  It  was  destroyed  by  sympa¬ 
thetic  detonation  during  cleanup  operations. 

SECTION  3 

ANALYSIS  OF  THE  GROUND  MOTION  DATA 


Up  to  now,  the  main  goals  of  data  analysis,  have  been  i)  to 
assess  the  credibility  of  the  data,  noting  where  ground  motion 
signals  were  affected  by  anomalies,  ii)  to  establish  gross  proper¬ 
ties  of  the  ground  motion  field  (velocity  of  first-arriving  signal, 
decay  of  peak  particle  velocity  with  range)  from  the  credible 
data,  and  iii),  by  examining  questionable  data,  to  determine  pro¬ 
bable  causes  of  errors  and  anomalies  along  with  means  of  removing 
them  -  thereby  recovering  useful  information. 

3.1  The  Data 

Data  from  the  ground-motion  instrumentation  array  was  recor¬ 
ded  in  analog  form.  Subsequently,  copies  of  the  master  analog 
tapes  were  used  to  generate  digital  records  covering  about  0.2  sec¬ 
ond  before  and  1.0  second  after  zero  time  for  each  event.  After 
analog-to-digital  conversion  and  recording,  plots  of  the  digital 
pulses  were  generated  (These  appear  in  Appendices  3  through  K  of 
Ref.  7). 

The  completed  ground  motion  instrumentation  array  could 
have  produced  680  velocity  gauge  records,  200  accelerometer  records 
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and  14  displacement  records  in  the  ten  events  of  Table  1.  Attri¬ 
tion  due  to  gauge  or  cable  failure,  connector  corrosion,  late 
installation  and  digital  record  processing  problems,  accounted  for 
the  loss  of  141  velocity,  68  acceleration  and  12  displacement 
records.  Pulses  representative  of  the  539  velocity  records,  132 
accelerations  and  4  displacement  gauge  records  reported  in  Ref.  7 
are  shown  in  Figures  3,  4  and  5.  The  pulses  in  Figure  3,  parts 
b-f,  show  that  noise  and  baseline  instability  were  significant 
problems  in  recording  velocities.  The  accelerometers,  Figure  4, 
show  significant  ringing  and  some  baseline  problems.  The  dis¬ 
placement  records,  Figure  5,  are  noisy  and  the  non-linear  cross¬ 
axis-coupled  response  of  the  gauge  adds  to  the  difficulty  of  their 
interpretation. 

Only  363  velocity  records,  87  accelerations  and  2  displace¬ 
ments  are  applicable  to  determining  the  characteristics  of  a 
spherically  symmetric  field  of  motion.  The  non-radial  velocity 
gauges  from  Event  III  can't  be  used  to  infer  radial  motion.  Like¬ 
wise,  since  signals  from  the  charge-pairs  cannot  be  separated 
reliably  into  single-charge  signals,  records  from  Events  IV-1.2, 
IV-2.2  and  IV-3  are  not  pertinent.  The  number  of  useful  records 
was  reduced  still  further  by  some  technical  problems  with  the  dig¬ 
ital  pulse  records.  These  problems  arose  during  attempts  to  read 
the  digital  ground  motion  records,  and  verify  them  by  comparison 
with  what  should  have  been  the  same  pulses  plotted  earlier  (Ref.  7) 
A  few  files  were  unreadable.  Identification  and  scaling  from 
"digital  bit  count"  to  measured  physical  amplitude  was  uncertain 
or  not  possible  for  many  others;  no  directory  was  provided  in  Ref¬ 
erence  7  for  the  tape  containing  most  of  the  digitized  accelera¬ 
tion  records.  This  left  the  plots  in  Ref.  7  as  the  primary  means 
of  defining  the  ground  motion  records.  Since  the  majority  of  the 
digitized  accelerometer  records  were  available  only  in  graphical 
form,  they  were  not  integrated  to  check  the  relative  accuracy  of 
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providing  signals  otherwise  good  through  peak  velocity,  could 
not  be  reliably  separated  into  single-charge  signals.  Determin¬ 
ation  of  peak  velocity  as  a  function  of  range  was  thus  limited 
to  the  117  good  records  (5  clean,  112  corrected)  from  single-charge 
events  listed  in  Table  3. 

Scale  factor  errors  were  present  in  plots  for  12  of  the 
good  records  for  single-charge  events  and  identification  ambigu¬ 
ities  affected  20  records  (one  of  them  also  affected  by  a  scale 
factor  error) .  The  identification  ambiguity  of  18  records  was 
resolved  by  reading  and  replotting  the  digital  records.  Scale 
factor  errors  consisting  of  transposition  of  digits,  substitution 
of  a  factor  for  a  different  gauge,  or  amplitude  shift  by  an  inte¬ 
gral  power  of  10,  were  corrected. 


3.1.2  Geometric  Correction  of  Data 

Since  the  velocity  gauge  array  was  designed  to  measure  the 
spherically  symmetric  field  of  motion  from  Event  III  and  latitudinal 
and  longitudinal  departures  from  symmetry,  the  motions  recorded 
for  all  other  events  must  be  adjusted  to  account  for  event-to-event 
changes  in  charge  location.  To  obtain  a  description  of  the  spher¬ 
ically  symmetric  radial  motion,  the  signal  measured  by  a  gauge 
pointing  in  some  arbitrary  direction  e^  must  be  related  to  the 

radial  velocity  V  =V  e  . 

r  r  r 

The  particle  velocity  at  any  point  is: 


V  =  V  e  +  V  e  +  V  e 
R  R  V  V  T  T 


(1) 


A  /\  A 

where  e  ,  e  ,  and  e  are  the  orthogonal  unit  vectors 
K  V  1 

the  directions  of  the  sensitive  axes  of  the  velocity 


pointing 
gauges . 
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The  spherical  radial  velocity  j_s 


vr  =  V  V  =  VvV  +  VvV  +  vT(VV 


where  e  is  a  unit  vector  pointing  along  the  ray  from  the  charge 
r 

to  the  gauge  canister.  However,  for  this  relationship  to  be  useful 
here,  records  for  all  three  components  of  motion  are  needed  and 
they  must  be  good  through  their  vector-sum's  peak.  Because  of 
uncertainties  and  omissions  in  the  data  linking  the  plots  of  Ref¬ 
erence  7  to  the  digitized  records  on  magnetic  tape,  the  construc¬ 
tion  of  vector  sums  was  not  attempted.  Instead,  it  was  assumed 


that  peaks  of  V_,VT_,  and  V  occur  simultaneously,  i.e.,  Eq.  (1) 
K  V  1 


is  assumed  valid  for  combining  component  peaks. 


Inverting  Equation  1,  leads  to  a  description  of  the  velocity 
component  measured  by  a  gauge  pointing  in  the  direction  e^  in 
terms  of  the  r,0,  and  ep-components  of  velocity: 


a  -* 
e.  .V 
i 


A  _  A  A  A  . 

e.*(V  e  +  V  e  +  V  e  ) 

r  r  9  9  p  cp 


=  Vvs>  +  VvV  +  VvV 


If  the  field  of  motion  is  spherically  symmetric  (V  =0  and 


V^=0),  Eq.  3  simplifies  to 


,A  A  . 

V .  =  V  (e.  . e  ) 
l  r  i  r 


Radial  velocity  can  then  be  obtained  as  follows  from  single  rec¬ 
ords  of  motion,  when  the  charge  and  gauge  location  are  known  along 
with  the  directions  of  the  sensitive  axes  of  the  gauges: 


V  =  V. /(e. .e  )  =  a.V. 
r  l  i  r  l 


/\  /\  —  2. 

where  *  =  (e^.e^)  is  the  geometric  correction  factor. 

If  the  field  of  motion  produced  by  a  single  charge  is 
spherically  symmetric  and  the  gauges  accurately  report  free-field 
motion,  then  all  the  estimates  of  the  radial  velocity  obtained 
from  a  given  canister  will  be  the  same.  For  a  canister  with  three 
gauges  (records  of  VR,VV,  and  VT),  Equation  2  will  yield  one  esti¬ 
mate  and  Equation  5  will  provide  three.  The  actual  field  will  be 
only  approximately  spherical,  the  gauges  only  approximately  aligned 
in  the  planned  directions,  the  canister  motion  a  bit  different 
than  that  of  the  free  field,  and  the  gauges  slightly  (and  differently) 
in  error  in  sensing  and  reporting  the  motion.  The  differences 
between  estimates  from  a  single  canister  will  reflect  the  combined 
effects  of  all  these  factors  (as  well  as  the  accuracy  of  our  - 
assumption  that  component  peaks  at  any  station  occur  simultaneously), 
and  will  give  a  qualitative  idea  of  measurement  accuracy.  Since 
several  canisters  contained  only  two  velocity  gauges  and  few 
canisters  with  3  gauges  yielded  complete  sets  of  good  records,  the 
estimates  of  radial  velocity  based  on  Equation  1  were  too  few  to 
be  of  value.  To  provide  an  adequate  data  base  and  allow  uniform 
treatment  of  the  velocity  data,  Equation  5  was  adopted  to  define 
the  amplitude  of  radial  velocity  for  each  single-charge  event. 

Because  the  description  of  canister  and  charge  locations 
given  in  Ref.  7  is  incomplete,  additional  data  were  sought  from 
records  of  telephone  conversations  and  working  notes  to  complete 
the  definition  of  gauge  and  source  geometry.  Although  some  of  the 
additional  data  led  to  conflicts  with  Ref.  7,  a  provisional  picture 
of  the  as-built  configuration  has  been  developed.  To  resolve  the 
conflicts  for  the  moment,  it  has  been  assumed  that  the  canister 
locations  of  Table  9  with  the  down-hole  survey  corrections  of 
Table  10  of  Ref.  7  are  correct  and  that  survey  layouts  for  the 
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slant  holes  '  provide  the  proper  orientation  of  the  slant  holes 
relative  to  the  rays  from  the  Event  III  charge  to  the  canisters. 
Where  data  are  missing**,  it  has  been  assumed  that  gauge  orienta¬ 
tion  and  location  are  correct,  i.e.,  unsurveyed  holes  did  not 
wander  from  their  intended  direction. 

Using  this  description  of  the  as-built  configuration,  the 
geometric  correction  factors  of  Eq.  5  were  computed  for  all  gauges 
for  all  events.  The  Event -III  factors  for  all  R-velocity  gauges 
were  compared  with  those  derived  from  gauge  mounting-block  details. 
This  comparison  revealed  aiming  errors  of  28°  for  gauges  D9UR  and 
D9UV,  1.8°  for  gauges  C7UR  and  C7UV,  and  an  indication  that  the 
inclination  of  in-plane  slant  hole  B9  was  7.2°.  Since  no  good 
records  were  obtained  from  B9  gauges,  no  further  effort  was  made  to 
resolve  the  question  of  slant-hole  inclination  or  gauge  orientation 
there.  The  1.8°  aiming  error  for  the  C7  gauges  was  considered 
small  enough  to  ignore  (errors  in  corrected  peak  radial  velocities 
from  this  source  would  be  less  than  7%,  and  less  than  47,  in  seven 
of  the  eight  cases  where  good  records  were  obtained).  Geometric 
correction  factors  for  gauges  D9UR  and  D9UV  were  recalculated  to 
account  for  their  misorientat ion. 

The  geometric  correction  factors,  estimates  of  radial  vel¬ 
ocity,  slant  ranges  and  yield-scaled  slant  ranges  are  listed  in 
Table  3.  The  peak  radial  velocity,  scaled  slant  range  data  for 
all  entries  with  geometric  correction  factors  (Eq.  5)  less  than 
5  (gauge-axis  pointing  within  78.46°  of  ray  from  charge)  have  been 
plotted  in  Figure  6.  Points  labelled  A  through  D  (A,  B  and  C  from 
Event  II-7  and  D  from  Event  IV-2.1)  are  clearly  outside  reasonable 
scatter  bounds  for  the  remainder  of  the  data.  They  have  been 
ignored  in  obtaining  the  7-event  regression  fit  shown. 

*The  two  sets  of  data  imply  azimuths  for  line  C,  Stations  3-11, 
that  differ  by  ~16°. 

v*Despite  diligent  search  and  a  careful  review  of  all  available 
data,  the  orientation  of  in-plane  slant  hole  B-9  remains 
unknown. 
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Figure  6.  Dependence  of  Peak  Velocity  on  Scaled  Slant  Range.  Data  from  97 
gauges  in  7  single-charge  events  gives  the  linear  regression  curve: 

UmaX=0*244R+  ‘ ^  ^  and  a  907,  confidence  band  whose  limits  are  in  the 

ratio  U  /U  =3.67.  Ignoring  the  less  consistent  data  from  Events  II-7  and 

IV-1.0  left  73  records  and  tne  5-Event  fit:  U  =  0.239R  1,^6±0.048  with 
,  max 

D  /U'=3.13. 


■Mf.v 


V  =  O^R-1-38340-050 
max 

The  90%  confidence  limits  cover  a  band  whose  bounds,  top  to  bottom, 
are  in  the  ratio  3.67:1.  Further  study  of  the  data  led  us  to  reject 
the  peaks  from  two  events  (Section  3.2  below),  giving  a  best  esti¬ 
mate  somewhat  different  than  this  one. 

3.2  Comparison  of  Single-Charge  Ground  Motion  Data  with  Cowboy 
and  Salmon. 

Ground  motions  were  measured  for  two  other  programs  con- 

2 

ducted  in  dome  salt:  Project  Cowboy  in  the  Carey  Salt  Mine  at 

Winnfield,  Louisiana,  and  Project  Dribble^-  in  the  Tatum  salt  dome 

near  Hattiesburg,  Mississippi.  Project  Cowboy  included  8  tamped 

shots  fired  between  17  December  1959  and  9  March  1960  that  were 

comparable  to  the  Cowboy  Trails  single-charge  events.  Explosive 

charges  in  the  series  of  tamped  shots  ranged  from  20  to  1000  lbs 
1 3 

of  Pelletol*  (spheroidal  pellets  of  TNT  with  diameters  of  about 
3/32  inch).  Under  the  conditions  for  tamped  charges  in  Cowboy  and 
Cowboy  Trails,  the  confinement -dependent  energy  release  of  Pelletol 
has  been  taken  as  3.4  kj/g,  the  value  deduced  from  the  earlier 

4 

phases  of  Cowboy  Trails. 


Data  for  events  4,  7,  9,  11,  16  and  17  with  charges  weighing 

100  to  1000  lbs  were  used  to  obtain  the  peak  velocity/range  curve 
a -1  527  14 

(Vmax=  0.1456R  ‘  )  shown  in  Figure  7.  Velocity-gauge  data  and 

integrated  accelerometer  data  from  the  Salmon  event  of  Project 

Dribble,  a  tamped  5.3  kt  nuclear  detonation  fired  22  October  1964 

a  _  I  q<524  14 

were  used  to  obtain  the  Salmon  curve  (V  =  0.1930R  '  ) 

max 


Careful  review  of  ground-motion  data  acquisition  systems 
used  in  Cowboy^  and  Salmon^-  has  shown  that  the  scaled  frequency 


*Trojan  Corporation's  Nitropel,  the  currently  available  pelletized 
TNT  explosive,  was  used  in  the  spherical  Cowboy  Trails  charges. 
DuPont  Pelletol  1  was  used  in  the  Cowboy  cylindrical  charges. 
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Figure  7.  Comparison  of  Linear  Regression  Fits  to  Peak  Velocity,  Range  Dat 
from  Cowboy  Trails,  Cowboy  and  Salmon.  Fil tered\Salmon  data  have  been  pass 
through  a  digital  filter  representing  the  W:i-scaled  frequency  response 

of  the  Cowboy  data  acquisi  ion  system.  High  frequency  response  is  reduced, 
decreasing  both  amplitude  and^the  rate  of  decay  with  range.  Curve  for  Cowb 
is  dashed  for  R>  1.61  km/kt  3  where  orientation  of  gauges  could  not  be 
confirmed. 


response  of  the  velocity  gauges  used  in  Cowboy  is  quite  low, 

especially  with  respect  to  that  for  Salmon.  Limited  frequency 

response  sharply  reduces  peaks  at  close-in  stations,  but  the  effect 

decreases  with  range;  higher  frequencies  decay  more  rapidly  with 

distance  than  lower  frequencies.  As  a  result  the  rate  of  decay 

of  neak  velocity  with  range  is  decreased  by  reducing  the  frequency 

response  of  the  instrumentation  system.  Passing  Salmon  data  through 

digital  filters  representing  the  frequency  response  of  the  Cowboy 

velocity  gauges  scaled  to  Salmon  yield  reduced  peak  velocities  by 

26  to  417,  close-in,  depending  on  the  gauge  type  represented  by  the 

16 

filter  and  by  6  to  237,  at  the  outermost  range .  The  decay-expon¬ 
ent  fell  from  1.9524  to  1.783  or  1.766  as  shown  by  the  curves 
bounding  the  shaded  region  for  filtered  Salmon  in  Figure  7.  No 
attempt  has  yet  been  made  to  account  for  the  differences  in  ground 
motion  from  cylindrical  (Cowboy)  and  spherical  (Cowboy  Trails) 
charges  or  for  the  differences  between  the  frequency  response  char¬ 
acteristics  of  velocity  gauges  used  to  record  Cowboy  and  Cowboy 
Trails  motions. 

The  exponent  n  for  decay  of  peak  velocity  with  range  (1.383) 
is  lower  in  Cowboy  Trails  than  that  observed-  for  both  Cowboy  (1.527) 
and  Salmon  (1.9524;  1.775±.009  after  filtering  to  compare  with 
Cowboy).  However,  907.  confidence  limits  encompass  the  regression 
curves  for  both  Salmon  and  Cowboy  over  the  scaled  range-interval 
from  .175  km/kt^  to  the  most  distant  Cowboy  gauge  (and  our  best 
present  estimate  of  n  is  somewhat  larger  than  1.383;  below).  Still, 
Cowboy  Trails’  peak  amplitudes  are  almost  all  higher  than  the  Cowboy 
regression  line,  which  lies  just  inside  the  lower  90%  confidence 
bound.  Thus,  while  little  overlap  occurs  in  the  scaled  range- 
intervals  on  which  ground  motion  was  observed  in  Salmon  and  Cowboy 
Trails  (0.37<Jl<0.42  km/kt3),  it  is  reassuring  to  find  that  the 
regression  line  for  Cowboy  Trails  falls  between  the  fits  to  Salmon 
and  filtered  Salmon  data. 


The  comparison  of  Cowboy  Trails  with  Cowboy  presents  some 
puzzles.  The  combination  of  higher  amplitude  but  lower  rate  of 
decay  with  range  cannot  be  explained  by  frequency  response  differ¬ 
ences  of  the  kind  demonstrated  in  comparisons  of  Cowboy  and  Salmon. 
The  differences  might  result  from  the  predominance  of  Cowboy  Trails 

observations  at  scaled  ranges  beyond  those  used  in  Cowboy  (the  log 

L 

mean  range  is  2.465  km/kt  ) ,  where  deformation  is  more  nearly  elas¬ 
tic,  or  from  systematic  errors  in  single  events.  In  particular, 
signals  from  Event  II-8,  the  shallowest  shot,  may  all  be  contamin¬ 
ated  by  surface  reflection;  records  for  Event  II-7 ,  the  deepest 
shot,  may  all  be  affected  by  a  suspected,  but  unconfirmed,  error 
in  charge  location  that  overstates  the  ranges  to  all  stations 
(Section  3.3).  The  full  set  of  records  may  also  be  affected  by 
superposed  noise  that  overstates  all  peak  values. 

Treating  good  records  from  each  event  separately  produces 
the  seven  regression  lines  shown  in  Figure  8,  whose  parameters  are 
listed  in  Table  4.  An  eighth  regression  line,  also  shown,  fits 
all  the  good  Cowboy  Trails  records  for  scaled  ranges  less  than 
1.61  km/kt3  -  the  largest  range  at  which  radial  orientation  of  the 
Cowboy  gauges  has  been  confirmed.  The  exponents  for  power-law 
decay  vary  from  a  low  of  1.232  for  Event  IV-1.0  to  a  high  of  1.775 
for  Event  II-7.  The  90%  confidence-limit  bards  have  widths,  U+/U  , 
varying  from  1.22  for  Event  III  to  4.63  for  Event  II-7.  The  decay 
exponents  for  Events  III,  IV-2.0,  IV-2.1  -  the  only  events  yield¬ 
ing  good  records  for  Rsl.7  km/kt ^  (the  interval  of  overlap  with 
Cowboy)  -  are,  respectively,  1.500,  1.522,  1.596  and  1.556.  For  all 
intents,  they  are  equal  to  that  for  Cowboy  (1.527).  Amplitudes, 
however,  run  from  about  1.20  to  nearly  twice  the  Cowboy  amplitude. 
Peak  velocities  from  the  four  events  yielding  records  only  for 
large  scaled  ranges  -  Events  II-?,  II-8,  IV-1.0,  and  IV“ 1*1  -  tend 
to  be  high.  Their  regression  lines  lie  above  those  for  Events  III, 
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Linear  Regression  Fit  Parameters  for 


••Omitted  from  5-Event  fit 


IV - 2 . 0 ,  and  IV-2.1. 


The  mass  of  data  from  these  events  (64  of  the  97  peak-velo¬ 
city/range  pairs)  biases  the  regression  line  toward  lower  decay 
rates.  The  regression  analysis  of  single  events  suggests  that  the 
decay  rate  over  the  larger  range  interval  may  be  as  high  as  that 
for  Cowboy.  It  also  shows  that  data  for  Events  II-7  and  IV-1.0 
are  less  consistent  than  data  for  the  other  events;  standard  devi¬ 
ations  are  two  to  four  times  the  next  largest  values.  On  the  basis 
of  those  large  uncertainties  and  the  increasing  evidence  that  the 
unsurveyed  charge  placement  hole  for  Event  II-7  wandered  far  from 
its  intended  location,  all  records  for  Events  II-7  and  IV-1.0  were 
ignored  and  a  revised  linear  regression  fit  derived.  The  revised 
fit  (shown  as  5-Event  Fit  in  Figure  6)  is 


U  =  0.2395R"1’  ±0,048° 

max 


The  standard  deviations  for  n  and  log  A  (the  decay  exponent  and 


logarithm  of  the  amplitude,  respectively)  and  the  width  of  the  90% 


confidence  bounds  (U  /U  )  are  reduced  from  the  seven-event  fit. 


To  test  the  Cowboy  Trails  data  for  transition  to  elastic 
decay  over  the  interval  of  observation,  quadratic  regression  fits 


log  U  =  A  +  B  log  R  +  C  (log  R)2 
°  max  vo/ 


were  obtained.  Results  for  Event  III  are  meaningless  since  only 
4  records  are  available  and  the  quadratic  fit  shows  decay  rate 
increasing  with  range  (C  >0) .  The  nine  other  cases  show  decay  rate 
decreasing  with  range,  but  the  uncertainties  in  A  and  B  [Eq.  (7)1 
(measured  by  the  standard  deviation  divided  by  the  coefficient) 
are  larger  than  the  corresponding  uncertainties  in  the  linear  fit. 
The  uncertainty  in  C  (standard  deviation  +  C)  ranges  from  0.48  to 
10.8.  The  large  standard  deviations  of  the  coefficients  show  only 


that  there  are  wide-ranging  quadratic  fits  of  almost  equal  merit 

(as  measured  by  the  least-squares  criterion).  To  get  a  better  iea 

of  the  significance  of  curvature,  the  standard  deviation  of  the 

narabola  as  an  estimator  of  U  was  compared  with  that  for  the 

max 

linear  fit  (corrected  for  the  change  in  the  number  of  degrees -of - 
freedom).  In  two  of  the  nine  relevant  cases  of  Table  4  (Event  III 
doesn't  count)  the  standard  deviation  was  greater  for  the  quadratic 
fit  than  the  linear  fit.  '  In  three  cases  it  was  smaller  by  less 
than  1%,  and  in  two  other  cases  by  only  1.970.  In  these  seven  in¬ 
stances,  the  fits  to  Events  II-7,  IV-1.1,  IV-1.0,  7-Events,  5-Events, 
and  Events  II-8,  and  IV-2.1,  respectively,  there  is  no  substantial 
evidence  that  the  quadratic  is  the  better  fit,  i.e.,  that  a  trans¬ 
ition  toward  elastic  decay  takes  place  as  range  increases.  In 
the  two  remaining  cases,  Event  IV-2.0  and  the  data  for  the  region 
of  overlap  with  Cowboy,  the  standard  deviation  for  the  quadratic 
showed  net  reductions  of  about  5  and  107o.  The  quadratic  fit  for 
the  Cowboy  overlap  region  has  an  equivalent  power  law  decay  (R  ) 
that  varies  from  n^.OS  at  the  inner  limit,  R=0.379  km/kt^,  to 
n=-0 . 01  at  R=1 . 583  km./kt^,  a  result  at  odds  with  elastic  decay 
(n=l)  beyond  the  range  interval  covered.  The  results  for  Event 
IV-2.0  (with  57o  smaller  standard  deviations  than  the  linear  fit) 

A 

are  consistent  with  transition  toward  elastic  decay  for  R  >  1.6 

jL  A 

km/kt  ,  the  equivalent  decay  rates  ranging  from  n=1.85  at  R=0.46 
1  ^  A. 

km/kt3  to  n=l . 15  at  R=8.47  km/kt3.  Overall,  however,  evidence 
for  a  decrease  in  n  with  slant  range  must  be  regarded  as  very  weak. 

Without  additional  data,-  linear  regression  descriptions 
of  the  decay  of  peak  velocity  with  range,  AR  ,  make  the 

most  sense.  Whether  the  data  are  treated  by  single  events  or 
groups  of  events,  the  decay  exponent  n  lies  between  the  value  ob¬ 
tained  from  Cowboy  data  (1.527)  and  about  1.35  (our  best  estimate 


being  1.45).  There  is  evidence  to  support  either  extreme  out  to  the 
largest  range  of  observation  here*  Also,  while  the  data  are  too 
scattered  to  determine  whether  the  decay  rate  decreases  with  range, 
decay  is  clearly  not  elastic  ( aR  even  at  scaled  ranges  as 
great  as  11.3  km/kt^,  or  seven  times  the  previous  maximum  scaled 
range  of  observation  from  tamped  shots. 

3.3  Propagation  Speed  of  First  Arrivals 

In  principle,  determining  first-arrival  times  is  simpler 
than  extracting  peak  material  velocities.  It  does  not  depend  on 
exact  knowledge  of  gauge  sensitivity  or  field  calibration  signal 
strength,  and  no  geometric  correction  factors  need  be  applied  to 
account  for  gauge  orientation.  All  that  is  required  is  i)  an  ' 
accurate  reference  timing  signal,  ii)  accurate  location  of  ex¬ 
plosive  charges  and  gauges,  and  iii)  an  identifiable  first  arrival 
signal.  Both  accelerometer  and  velocity  gauge  signals  can  be  used 
and  the  gauges  need  not  function  flawlessly  after  first  arrival 
(till  peak  velocity  or  later).  In  addition,  both  single-charge 
and  charge-pair  events  can  be  used. 

In  practice,  noise  often  made  first-arrival  detection  dif¬ 
ficult  in  the  Cowboy  Trails  data.  More  frustrating  was  the  occas¬ 
ional  loss  of  synchronization  with  the  simultaneously  recorded 
IRIGB  timing  signal  which  introduced  unknown  time  shifts  into  a 
number  of  record  plots.  In  spite  of  these  difficulties,  times - 
of-arrival  were  determined  for  368  Cowboy  Trails  records,  including 
records  from  all  events  (see  Table  5).  In  several  cases,  the  data 
suggest  otherwise-undetected  time-shifts  and/or  appreciable 
uncertainties  in  gauge  and  charge-location.  To  exclude  the  most 
obvious  of  these,  a  wavespeed  C'  was  estimated  from 

C*  =  r/(t  -T  ) 
a  o 
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If  this  estimate  was  between  3C  / 4  and  4C  / 3 ,  where  C  =  4.6  m/ms 

oo  o 

(the  expected  wave  speed),  the  data  were  retained. 

Difficulties  encountered  in  drilling  the  charge  emplacement 
hole  and  reaming  the  spherical  charge  cavity  for  Event  II-7  sugges¬ 
ted  that  the  hole  had  wandered  significantly  from  vertical.  Dig¬ 
ital  data  file  identification  problems  had  cast  doubt  on  the  lo¬ 
cations  (through  gauge  identification)  of  most  of  the  records  for 
Evept  IV-3  and  a  few  isolated  records  in  other  events.  As  a 
further  test  of  the  validity  of  the  data,  regression  fits  to  range 

as  a  function  of  arrival  time  r=r  +  C(t  -T  )  were  obtained  for 

o  a  o 

each  event,  Table  5.  Large  values  of  rQ  suggest  systematic  errors 
in  reported  values  of  the  slant  range  or  time  of  arrival.  Large 
values  of  the  standard  deviations  of  r,  rQ  and  C  suggest  random 
errors.  Data  from  Event  II-7  is  suspect  on  both  grounds.  It  also 
shows  a  suspiciously  low  wave  speed  4.07  m/ms.  All  Event  II-7 
time-of-arrival  data  was  therefore  excluded  from  further  analysis. 

Data  from  Event  IV-1.2  also  shows  larger  standard  deviations 
and  a  larger-than-average  .  Since  there  was  no  other  reason  to 
be  wary  of  the  data  for  this  event,  it  was  retained.  The  resulting 
data  set  -  292  time-of-arrival,  range  pairs  from  nine  events  -  was 
used  to  determine  the  propagation  speed  for  first  arrivals,  4.479 
±0.021  m/ms.  The  regression  fit  and  all  368  data  points  are 
shown  in  Figure  9. 

SECTION  4 

ANOMALOUS  VELOCITY  GAUGE  SIGNALS 

The  completed  Cowboy  Trails  instrumentation  array  should 
have  provided  434  velocity  gauge  records  suitable  for  determining 


TIME  OF  ARRIVAL,  ms 

Figure  9.  Cowboy  Trails:  Time-of -Arrival ,  Range  Data.  Ihe  linear  regression 
line,  r  =  r0  +  C0ta,  is  the  least-squares  fit  to  292  points  from  9  of  the  10 
events.  It  indicates  a  first-arrival  signal  speed  of  4.479  ±  0.021  m/ras. 
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radial  velocity  from  the  seven  single-charge  events 
of  those  records  (347)  exhibited  noise  and  baseline 
Table  2),  56%  of  them  being  so  severely  affected  th 
considered  valid  through  peak  velocity.  Outright  f 
gauges  accounted  for  the  loss  of  only  75  records, 
number  lost  to  noise  and  baseline  problems  is  more 
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Figure  10.  Simplified  Schematic  Diagram  of  the  Bell  &  Howell 
Velocity  Gauge.  A  pulse  generator  was  substituted  for  the 
piezoelectric  crystal  in  most  of  the  bench  testing. 
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to  Recording  Van 


Figure  11.  Schematic  of  the  Bias  Nulling  Circuit  for  the 
Bell  &  Howell  Velocity  Gauge. 


circuit  was  driven  by  a  pulse  generator  in  place  of  the  piezoelec¬ 
tric  accelerometer  element  (crystal)  of  the  gauge.  Once  a  leaky 
capacitor  had  been  identified  and  replaced,  the  breadboard  circuit 
performed  as  predicted  from  the  circuit  equations. 

Laboratory  testing  of  the  unpotted  gauge  required  the  addi¬ 
tion  of  the  bias  nulling  circuit  (Fig.  11)  used  in  the  field  to 
remove  a  12  volt  DC  signal  superimposed  on  the  dynamic  output  of 
the  integrator  circuit.  The  pulse  generator  was  also  used  to 
drive  the  factory  integrator  circuit.  When  the  gain  of  the  nulling 
circuit  was  accounted  for,  the  factory  integrating  circuit  per¬ 
formed  according  to  specifications,  matching  (within  a  few  percent) 
the  breadboard  response  and  calculated  response  of  the  circuit 
equations  to  steady  state  sinusoidal  input  of  from  1  to  5000  HZ 
as  shown  in  Table  6.  Higher  gains  but  comparable  agreement  with 
calculated  response  were  obtained  for  single  boxcar  pulse  input, 
for  pulse  durations  from  1  ms  to  a  few  hundred  ms . 

The  integrator  circuit  was  also  tested  by  attaching  the 
piezoelectric  accelerometer,  shock  loading  it  while  monitoring 
both  the  accelerometer  and  integrator  circuit  output.  The 
acceleromecer  output  pulse  (Fig.  12a)  was  then  modelled  by 
superposing  19  boxcar  pulses  and  used  to  drive  the  circuit 
equations.  Calculated  and  measured  integrator  circuit  output 
pulses  are  compared  in  Figure  12b. 

The  integrator  circuit  response  was  a  linear  function  of 
input  pulse  amplitude  up  to  the  limiting  voltage  specified  by  the 
manufacturer.  At  higher  input  voltages  (accelerations  beyond  the 
certified  operating  limit)  anomalous  integrator  output  was  seen. 

The  anomalies,  however,  were  very  different  from  the  behavior 
noted  in  the  Cowboy  Trails  data.  The  accelerometer  output  voltage 
and  bias -nulled  integrator  output  voltages  for  an  overload  test 
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*These  are  single  measurements.  Other  entries  for  breadboard  and  factory 
circuits  are  averages  of  two  measurements, 


a)  Accelerometer  Stage  Output  (Integrator  Input) 


✓  i 

!  A 


b)  Velocity  Pulses  (Integrator  Output) 


Figure  12.  Comparison  of  Observed  and  Calculated  response  of  a 
Bell  &  Howell  Velocity  Gauge  to  an  Observed  Acceleration  Pulse. 


are  shown  in  Figure  13.  When  the  accelerometer  crystal  output 
voltage  first  exceeds  10  volts  (the  saturation  voltage  for  the 
unity  gain  amplifier  between  the  crystal  and  the  integrator  cir¬ 
cuit  in  Figure  10)  the  integrator  circuit  abruptly  ceases  to 
function  properly,  point  A.  The  subsequent  abrupt  flattening  of 
the  accelerometer  output  at  ~12  volts  occurs  too  late  to  have 
triggered  the  malfunction.  It  has  no  visible  effect  on  integrator 
response.  The  integrator  circuit  returns  to  proper  function 
when  the  accelerometer  voltage  drops  below  10  volts,  point  B,  but 
has  lost  its  reference.  The  output  signal  is  meaningless  beyond 
point  A. 

The  accelerometer  stage  was  also  subjected  to  cross -axis 
loading  -  an  apparently  reliable  means  of  exciting  baseline  drift 
in  NMERl's  evaluation  tests'^.  A  small  spring-driven  jiggle 
table  carrying  the  piezoelectric  crystal  was  moved  to  one  of  its 
extreme  positions  and  released,  subjecting  the  accelerometer  stage 
to  a  slowly  decaying  ~30  HZ  signal  with  a  peak  acceleration  of 
1  or  2  g's,  much  lower  than  the  nominal  200  to  1000  g’s  of  the 
NMERI  tests.  Velocity  signals  for  both  onraxis  and  cross-axis 
loading  wer  clean,  sinusoidal  pulses  without  baseline  drift. 
Amplitude  of  the  output  for  cross -axis  loading  was  about  10%  of 
that  for  on-axis  loading. 

Based  on  the  long  series  of  bench  tests  on  factory  and  bread 
board  inttgrator  circuits,  the  integrator  is  clearly  not  the  cause 
of  baseline  problems.  Within  the  specified  operating  limits,  its 
steady  state  response  to  sinusoidal  input  is  that  computed  for 
the  circuit  (phase  and  gain).  The  circuit  is  stable,  linear,  and 
accurate  for  single  boxcar-pulse  and  pulse-train  input  for  dura¬ 
tions  from  0.5  to  a  few  hundred  milliseconds  and  input  signal  ampli 
tudes  less  than  10  volts.  If  baseline  drift  is  the  result  of  a 


gauge  defect,  the  defect  is  not  in  the  integrator  circuit.  Tests 
of  the  accelerometer  stage  were  not  complete  enough  to  rule  it 
out  as  the  source  of  baseline  drift  problems;  amplitudes  were  much 
too  low.  Exposure  to  the  dank,  salty  atmosphere  of  the  mine  which 
corrosively  attacked  tools,  electronic  gear  and  unprotected  circuit 
connections,  may  have  contributed  to  the  baseline  problems  of  the 
velocity  gauge*,  but  the  results  of  the  NMERI  tests  (not  explained) 
in  the  benign  environment  of  the  laboratory  suggest  otherwise. 

4.2  Noise 

During  screening  and  evaluation  of  the  ground  motion  data 
for  Cowboy  Trails,  noise  frequencies  were  estimated  by  counting 
the  number  of  peaks  over  time  intervals  of  a  few  tens  of  milli¬ 
seconds.  Since  the  noise  signals  are  not  clean  sinusoids,  the 
frequencies  are  only  approximate.  High  voltage  spikes  in  the  com¬ 
mercial  60  HZ  AC  power  (the  result  of  a  load  switching  surge) 
triggered  the  early  firing  of  the  IV-2a  charge  (Event  IV-2.0).  As 
a  result,  records  were  expected  to  show  60  HZ  noise  and  its  higher 
harmonics .  Spikey  60  HZ  noise  (see  Figure  14)  that  nearly  obscures 
the  ground  motion  signal  was  found  in  some  records. 

Examination  of  the  records  show  noise  frequencies  of  300 
to  600  HZ  as  well  as  the  expected  60  HZ  already  noted.  Some 
records  show  a  ringing  response  to  signal  arrival,  suggesting 
vibration  of  a  gauge  mount  or  gauge  element.  Noise  bursts  and 
ringing  signals  showed  frequencies  of  about  1500  HZ.  In  those 
instances  where  a  regular  pattern  of  noise  can  be  identified,  it 
appears  that  signal  quality  can  be  enhanced  by  further  processing 
of  the  digital  records  .  .  e.g.,  the  use  of  notch  filters  to  remove 
narrow-band  noise.  More  complex  procedures  (signal  differencing  in 


Figure  14.  Spikey  60  HZ  Noise 
this  accelerogram,  the  noise  S] 
pulse  which  arrives  at  about  3 


transform  space)  may  be  required  to  remove  the  spikey  60  HZ  noise 
pulses  like  those  seen  in  Fig.  14. 

The  techniques  for  removing  noise  don't  require  that  the 
sources  of  the  noise  be  identified,  though  that  would  be  helpful 
in  preventing  data  loss  to  noise  in  future  events.  The  suspected 
sources  of  noise  here  are  the  commerical  AC  power,  shock-induced 
vibration  of  gauge  mount  and  gauge  elements,  corrosion,  and  diff¬ 
erences  in  the  ground  potentials  in  the  mine  and  at  the  instru¬ 
mentation  van  at  the  ground  surface. 


SECTION  5 

CONCLUSIONS  AND  RECOMMENDATIONS 


Ground  motion  data  from  single  charge  events  in  the  Cowboy 
Trails  program  show  that  signal  propagation  in  dome  salt  is  inelas 
tic  out  to  scaled  ranges  of  11.3  km/kt^,  a  range  seven  times  that 
where  good  free-field  measurement  had  been  made  before.  Over  the 

/S  i 

scaled  range  interval  0.388  <  R<  11.3  km/kt  peak  velocity  decayed 
according  to 


V 

max 


0.2395R 


-1.456±0. 048  . 

m/s 


The  rate  of  decay  is  slightly  lower  than  that  observed  in  Cowboy, 
but  amplitudes  are  higher  by  about  a  factor  of  two. 

Propagation  speed  for  first  arrivals  from  both  single-charge 

A 

and  charge-pair  events  over  the  scaled  range  interval  0.388  <  R 


<  9.76  km/kt3  was 


Cq=  4.479*0.021  km/s 

in  close  agreement  with  the  values  for  Cowboy  (4.5  km/s  average^ 
from  shot  data,  4. 374*. 03  km/s  from  seismic  survey^)  and  Salmon 
(4 . 699*. 025  km/s  from  shot  data,  4. 507*. 135  from  uphole  survey, 
and  4.516*.  045  from  cross-hole  survey^"). 

These  results  are  based  on  a  small  fraction  of  the  records 
obtained  .  .  73  of  514  velocity  records  and  292  of  625  times-of- 
arrival.  Analysis  of  subsets  of  data  (i.e.,  data  for  separate 
events)  led  to  the  exclusion  of  all  time-of-arrival  and  velocity 
data  for  Event  II-7 .  A  large  error  in  source  location,  hinted  at 
by  difficulties  in  the  drilling  and  reaming  operations,  was  appar¬ 
ent  in  the  time-of-arrival  analysis  and  in  the  fit  to  peak  velocity, 
range  data.  Peak  velocity  data  from  Event  IV-1.0,  the  result  of 
the  first  misfire  of  a  charge-pair  event,  was  also  excluded  on  the 
basis  of  standard  deviations  substantially  larger  than  those  for 
fits  to  other  events.  Baseline  drift  and  noise  were  severe 
enough  in  325  records  to  render  them  invalid  for  determining  peak 
velocity.  The  inability  to  reliably  separate  pulses  from  the  two 
charges  in  Event  IV-1.2,  IV-2.2  and  TV-3  and  large  apparent  errors 
in  the  locations  of  charges  IV-2.2a  and  IV-3a  led  us  to  forgo  use 
of  charge-pair  data  in  assessing  peak  velocities  from  single  charges 
Ambiguity  in  the  azimuthal  coordinate  for  canisters  C3  through  Cll 
clouded  the  data  from  those  gauges  in  all  but  Event  III  through  an 
orientation  uncertainty  of  about  16  degrees.  Locations  for  those 
canisters  were  uncertain  by  about  20%  of  slant  range  from  Event  III 
charge.  In  spite  of  all  these  difficulties  the  data  show  a  high 
degree  of  consistency  in  both  peak  velocity  as  a  function  of  range 


and  first-arrival  speed.  It  therefore  seems  clear  that  further 
effort  to  recover  data  from  the  flawed  records,  and  to  reduce 
uncertainties  in  configuration,  would  be  worthwhile. 


Specific  recommendations  follow: 

1.  Perform  laboratory  tests  of  the  piezoelectric  integrating  accel¬ 
erometer  to  determine  the  cause  of  baseline  drift. 

2.  Correct  velocity  pulses  from  the  300  records  affected  by  base¬ 
line  problems,  using  results  of  the  baseline  drift  tests  in  con¬ 
junction  with  digital  data  processing  techniques. 

3.  Try  removing  noise  from  the  accelerometer  and  velocity  gauge 
records  by  such  means  as  spectral  analysis  and  narrow-band  filtering. 

4.  Reduce  the  uncertainty  in  the  location  of  the  charge  for  Event 
II-7  inverting  time-of-arrival  data,  using  the  first  arrival  wave- 
speed  and  gauge  stations  with  well-established  positions.  The  slant - 
range  uncertainty  for  canisters  along  C-line  is  probably  too  small 

(<  10  m)  to  be  resolved  in  a  similar  way. 

5.  Integrate  the  accelerometer  records  correcting  first  for  non- 
spherical  orientation.  Compare  with  corresponding  velocity  records 
to  establish  the  relative  performance  of  the  two  gauges  (acceler¬ 
ometer  and  integrating  accelerometer)  if  unambiguous  access  to  the 
digital  records  can  be  obtained. 

6.  Redigitize  velocity  and  accelerometer  data  in  events  where 
many  records  remain  inaccessible  due  to  uncertain  identification, 
uncertain  reference  time,  or  digital  tape  read-write  problems. 

7.  When  corrections  for  baseline,  noise  and  as -built  geo¬ 
metry  have  been  completed  and  the  technical  problems  of  digital 
record  identification  and  access  resolved,  pass  the  Cowboy  Trails 
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records  through  filters  representing  the  Cowboy  gauge  frequency 
response  characteristics. 

8.  Continue  to  develop  displacement  gauges.  Accurate  measure¬ 
ment  of  displacement  at  one  or  more  times  establishes  a  bench¬ 
mark  for  correcting  baseline  problems  in  velocity  and  acceleration 
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APPENDIX  A 


THE  PLANNED  COWBOY  TRAILS  PROGRAM 

A.l  Program  Objectives 

The  program's  primary  objectives  were  to 

a)  Measure  a  nearly-spherical  seismic  source  driven 
by  a  contained  explosion; 

b)  Determine  an  actual  spectrum  of  NE-CE  equivalence 
factors ; 

c)  Determine  the  site-to-site  variability  of  explosi¬ 
vely-driven  motion  in  "one"  medium; 

d)  Determine  the  accuracy  of  simple  scaling  rules  from 
lab  yields  (10  ^*^kt)  to  that  of  Salmon  (10'^kt). 

The  program  also  had  a  set  of  secondary  objectives,  namely,  to 

e)  Measure  the  variation  of  the  field  with  overburden; 

f)  Measure  the  energies  of  common  and  promising  CE's; 

g)  Measure  the  effect  of  charge-shape  on  the  field; 

h)  Measure  permanent  displacement,  Dm,  vs.  slant  range; 
from  those  data,  determine  the  actual  drift  of  ground- 
motion  gauges  and  show  how  to  correct  for  it. 

In  addition,  the  program  served  almost  inevitably  some  ancil¬ 
lary  objectives  that  hold  broad  interest;  these  were  to 

i)  Test  low-cost  CE's  for  reproducibility; 

j)  Determine  how  feasible  it  is  to  use  cavities  formed 
by  CE  in  place  of  NE-formed  cavities,  as  charge-holes; 

k)  Determine  the  suitability  of  explosively-formed 
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ities  as  charge-holes  for  tamped  and  decoupled  shots 
("cheap  mining") ; 

1)  Test  stress  gauges. 

A. 2  The  Program 

The  program  visualized  to  gain  these  objectives,  which  are  all 
experimental ,  hinges  on  a  series  of  experiments  requiring  only  post 
shot  measurements  (Phases  I  and  II  of  the  program)  and  another  ser¬ 
ies  in  which  motion  is  measured  dynamically  (Phases  III  and  IV). 

I.  The  following  scheme  was  proposed  to  determine  the  influence 

of  major  variables  on  cavity  size: 

a.  Charge  Shape  -  Compare  volumes  of  cavities  from  lab  shots 
with  cylindrical  charges  of  average  COWBOY  shape,  and  with  spher¬ 
ical  charges,  in  pressed  salt  and  also  in  blocks  of  Grand  Saline 
salt.  If  the  need  arises,  fire  a  COWBOY-like  charge  weighing 
~180  lbs  in  the  Grand  Saline  dome,  as  an  addition  to  Phase  II. 

b.  Overburden  -  Compare  the  shapes  and  sizes  of  Phase  II  cav¬ 
ities;  overburden  is  to  vary  from  40  to  100  bars  in  Phase  II.  Also, 
compare  the  shapes  and  sizes  of  cavities  from  lab  shots  in  pressed 
salt,  and  again  from  lab  shots  in  blocks  of  Grand  Saline  salt  (60- 
to-350  bars  of  confining  pressure). 

c.  Site  -  Compare  COWBOY  and  Phase  II  cavities,  corrected  for 
charge-shape  (la,  above),  to  get  a  direct  measure  of  the  differences 
between  salt  domes  as  explosion  sites  (overburden,  yield  and  even 
explosive ,  can  be  held  constant  in  that  comparison).  Compare  cav¬ 
ities  left  by  lab  shots  in  pressed  salt,  with  those  left  by  lab 
®hots  in  blocks  of  Grand  Saline  salt.  Compare  the  latter  cavities 


with  Phase  I  cavities  to  detect  differences  between  laboratory 
and  in-situ  behavior  of  Grand  Saline  salt. 


d.  Scale  -  Compare  cavities  from  Phases  I  and  II,  correcting 
(if  necessary)  for  the  small  difference  in  overburden  (lb,  above). 

II.  Influence  of  Major  Variables  on  Free-Field  Motion: 

a.  Charge  Shape  -  In  the  laboratory,  fire  cylindrical  charges 
of  average  COWBOY  shape,  and  spherical  charges,  in  both  pressed 
salt  and  blocks  of  Grand  Saline  salt.  If  the  need  arises,  fire  a 
charge  of  each  shape  in  its  own  block  of  dome  salt  -  or  even  a 
COWBOY-like  cylindrical  charge  (~180  lbs),  at  Grand  Saline. 

b.  Overburden  -  Compare  motions  from  Phase  III  and  a  deep 
shot  in  Phase  II;  the  overburden  variation  should  be  ~40  bars. 
Compare  motions  from  lab  shots  in  pressed  salt  at  confining  press¬ 
ures  from  0  to  50  bars,  and  again  at  confining  pressures  of  40,  100 
and  180  bars  in  blocks  of  Grand  Saline  salt. 

c.  Site  -  Compare  COWBOY  and  Phase  III  motions,  corrected  for 
charge  shape  (Ila,  above),  to  obtain  a  direct  comparison  between 
salt  domes  (yield  and  overburden  are  very  nearly  the  same  for 
COWBOY  and  Phase  III).  Also,  fire  lab  shots  in  pressed  salt  and 
dome-salt,  to  see  how  motions  differ  in  salt  media  created  by 
different  processes. 

d.  Material  Hardening  -  Compare  motions  measured  in  Phases 

III  and  IV  to  detect  any  systematic  change  in  the  free  field  driven 
by  a  single  spherical  charge,  as  the  number  of  previous  shots  grows 
(e.g.,  the  medium  could  become  noticeably  more  elastic). 

e.  Scale  -  Lab  shots  in  blocks  of  Grand  Saline  salt  (6-10  gms 
of  CE) ,  the  Phase  III  shot  in  the  Grand  Saline  dome  (180  lbs  of 
CE),  and  SALMON  (5.3  kt;  NE)  determine  the  effects  of  scale  on  free 
field  motion  over  nine  orders  of  magnitude  in  yield.  The  lab  and 
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Phase  III  shots  are  fired  in  the  same  medium  at  the  same  over¬ 
burden  (~60  bars),  but  the  SALMON  shot  was  fired  in  the  Tatum  dome 
at  an  overburden  of  ~180  bars ;  corrections  for  those  differences 
are  made  to  SALMON  data  using  results  obtained  in  this  program  on 
the  effects  of  site  and  overburden  variations  (lib.  and  c.,  above). 
Little  would  be  gained  by  comparing  adjusted  SALMON  pulses  in  all 
detail  with  lab  pulses  and  Phase  III  pulses  ,  because  NE  and  CE  dif¬ 
fer  too  greatly.  Instead,  in  comparing  fields  driven  by  different 
types  of  explosive,  define  scale  effects  in  terms  (for  example)  of 
rates  of  decay  of  peak  velocity  and  displacement  with  respect  to 
slant-range.  SALMON  and  COWBOY  pulse-shapes  proved  similar  enough  to 
justify  such  comparisons. 

III.  Seismic  Sources  and  NE-CE  Equivalence. 

a.  Seismic  Source  Created  by  NE  -  From  the  pulses  of  Phase 
III,  construct  a  complete  spherical  velocity  field  for  that  event 
over  its  instrumented  slant -range-interval  (the  "baseline  spheri¬ 
cal  field").  Do  likewise  for  SALMON,  and  extrapolate  the  SALMON 
field  outward  in  range  until  the  lowest  Phase  III  amplitudes  are 
reached;  to  extrapolate,  apply  to  SALMON  the  factors  that  define 
changes  in  pulse-shape  and  amplitude  in  the  baseline  field,  during 
whose  construction  those  factors  are  found.  Then,  after  correcting 
superposition-data  from  Phase  IV  for  the  effects  of  prior  shots, 
identify  a  minimum  slant -range,  or  "linear  radius",  beyond  which 
waves  from  a  single  CE  charge  propagate  in  linear  fashion.  As 
the  linear  radius  for  SALMON,  take  the  slant-range  at  which  both  (i) 
the  amplitude  of  the  SALMON  velocity  pulses  is  at  least  as  low  as  in 
the  Phase  III  field  at  its  linear  radius,  and  (ii)  over  the  fre¬ 
quency  band  of  interest  in  nuclear  monitoring,  logarithms  of  RVP- 
spectral  amplitudes  decay  at  the  same  rates  in  the  two  fields. 
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b.  Seismic  Source  Created  by  CE  -  Phase  III  motion  at  the 
linear  radius  found  in  Phase  IV  (Ilia,  above)  is.  the  seismic 
source  for  chemical  explosions  in  salt. 

c.  NE-CE  Equivalence  -  Bring  the  baseline  spherical  field 
(Ilia,  above)  and  the  extended  SALMON  field  to  a  common  yield, 
overburden  and  dome  by  adjusting  either  field  or  both  for  changes 
in  scale,  overburden  and  site  (the  adjustment  is  made  in  accord 
with  previous  results;  lib ,  c  and  e,  above).  From  RVP-spectra  of 
the  resulting  CE-  and  NE-driven  fields,  compute  equivalence-factor 
spectra  for  both  phase  and  amplitude.  Do  likewise  for  other  yields, 
burial  depths  and  domes,  if  the  effects  of  changes  in  scale,  over¬ 
burden  or  site  on  the  simply-scaled  baseline  field  have  been 

found  significant. 


A. 3 

A  flow  diagram  showing  how  the  measurements  made  in  each  phase 
of  the  program  contribute  to  achieving  its  goals  is  shown  below 
in  Fig.  A.  To  explain  how  the  diagram  works,  suppose  one  asks  how 
the  effects  of  scale  on  the  free  field  are  determined  in  the  pro¬ 
gram.  The  answer  is  found  by  tracing  backward  from  the  input  to 
"SCALE"  (for  the  free  field),  which  quickly  branches  in  three  di¬ 
rections.  One  branch  leads  back  to  "BASELINE  SPHERICAL  FIELD" 
[obtained  in  Phase  III]  and  "PELLETOL  ENERGY"  [obtained  by  combin¬ 
ing  data  from  Phases  I  and  II].  Another  branch  leads  directly  back 
to  the  FREE  FIELD  obtained  for  DOME  SALT  in  Lab  Shots.  The  third 
branch  leads  back  to  a  node;  there  data  from  SALMON  are  combined 
with  data  on  the  change  in  free  field  with  SITE  and  OVERBURDEN. 
Thus,  shots  at  lab  scale  (6-10  gms  of  CE) ,  mine  scale  (180  lbs  of 
CE)  and  the  Salmon  shot  (5.3  kt;  NE)  are  the  sources  of  data  as 
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